Persistent currents in multicomponent Tomonaga-Luttinger liquid: application to 
mesoscopic semiconductor ring with spin-orbit interaction 
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We study persistent currents in semiconductor ballistic rings with spin-orbit Rashba interaction. 
We use as a working model the multicomponent Tomonaga-Luttinger liquid which arises due to 
the nonparabolic dispersion relations of electrons in the rings with rather strong spin-orbit cou- 
pling. This approach predicts some new characteristic features of persistent currents, which may be 
observed in experimental studies of semiconductor ballistic rings. 
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I. INTRODUCTION 

It has been known for a long time that due to the 
conservation of the electron phase coherence an isolated 
metallic mesoscopic ring threaded by a magnetic flux may 
carry a persistent current^. This quantity, which equals 
to the derivative of free energy with respect to the flux, 
has been theoretically shown to exist in both ballistic 
ringa^ and rings with disordcr'^j. The effect of impuri- 
ties on persistent currents has been studied later on in a 
number of theoretical work a'^i^i^i'^i^i^ . 

The experimentally measured currentsiSiii*i^ in dif- 
fusive metallic rings demonstrated a disagreement with 
the theoretically predicted values. This observation 
has initiated an intensive debate. Several authors have 
studied the effect of Coulomb interaction on persistent 
current^i^'^i-^'^i^^i^^i^'^'i-^^i^^-^", either without or with dis- 
order. There is still no total consensus on the simul- 
taneous effect of impurities and electron-electron (e-e) 
interactions on persistent current. 

On the other hand, the experimental results for bal- 
listic GaAs/GaAlAs-based rings^i were found to be in a 
good agreement with the theoretical predictions. 

For one-dimensional (ID) systems there is a variety of 
powerful methods which allow to treat e-e interactions. 
For instance, the persistent current in the Hubbard ring 
with a flux and Coulomb repulsion can be found from the 
Bethe Ansatz (BA) solutioni^S^. The effects of disorder 
in the Hubbard model have been studied in Ref. [I^ by 
means of the bosonization technique. 

It is generally believed that the low-energy physics of 
ID systems can be often described using the Tomonaga- 
Luttinger (TL) Hquid concepts^ (see Refs. |24,25 for the 
review). In particular, quantum wires made of semicon- 
ductor heterostructures have been found to demonstrate 
the TL features^^. Therefore, the TL model appears to 
be a reasonable approximation for accounting for differ- 
ent effects in the ID semiconductor heterostructures. 

In Ref. (23 the persistent current has been studied 
within the TL model, and the particular emphasis has 
been made on the effects of electron number parity. The 
crucial significance of parity effects for the persistent cur- 



rent in ID systems has been recognized quite a long time 
agoiii. In particular, it has been noticed^^ that for a sys- 
tem of Ne interacting spinless fermions in the TL phase 
the sign of a current, which reflects the properties of the 
ground state, is diamagnetic for odd iVg, while it is para- 
magnetic for even Nf.. 

For the case of spinful interacting fermions the na- 
ture of the ground state may change depending both on 
the strength of Coulomb repulsion and on the relative 
parity of spin-up and spin-down electron numbers. For 
instancei^, in the Hubbard rings with an odd number of 
electrons N^. = N^.] + A^ei the current is paramagnetic 
around $ = 0, and its period is half the flux quantum 
$0 = hc/e. Meanwhile, in the Hubbard rings with even 
A"e and low filling the current is diamagnetic near zero 
flux, whereas at densities close to half-filling the currents 
may become paramagnetic. 

In the modern expanding field of spintronics the main 
research interest focuses basically on spin injection and 
spin detection in solid-state device922i2£. Semiconductor 
InAlAs/InGaAs-based heterostructures provide promis- 
ing opportunities for spin manipulation, because in such 

30,31 



materials spin-orbit (SO) splitting of Rashba type^i^'^'i is 
rather large, and this allows for especially strong coupling 
of spin polarization and electric field. A naturally aris- 
ing problem is to understand the effect of SO coupling 
on such a ground-state property as persistent current in 
a mesoscopic ring made of semiconductor material. In- 
deed, the issue of SO effects in noninteract ing, ID ba llistic 
rings has been already addressed in Refs. 32 33 34. Sur- 
prisingly, the exact BA solution which takes into account 
the SO effects in the appropriately modified Hubbard 
model is also available^^. Persistent currents in this sys- 
tem have been analyzed on the ground of the B A solution 
in Ref. l3a . The interplay of Coulomb repulsion and the 
SO coupling in combination with the parity effects leads 
to the remarkable features in the nature of current: al- 
though the SO interaction does not produce any effect in 
the case of a ring with the odd number N,^, the result for 
even-Ae rings differs from that of the standard Hubbard 
model without SO interaction. 

In the regime of low electron densities the ID sys- 



tern of interacting electrons with SO coupling can also 
be described in terms of a multicomponent TL model. 
The persistent current in the SO-splitted TL liguid has 
been studied in Ref. |33 using the approach developed in 
Refs. i27ii3& 

We have been discussing so far the effects of SO cou- 
pling which are caused by the relative horizontal shift 
of spin-up and spin-down parabolic dispersion curves. In 
such a situation the densities of states near the Fermi en- 
ergy are the same in each branch, and therefore they can 
be characterized by a unique value of the Fermi velocity. 
However, it has been recently shown^Si^SiiS that the pre- 
sence of strong SO interaction can qualitatively change 
the band structure of a 2D electron gas confined to the 
ID geometry. This can be achieved by applying both a 
specific lateral confining potential in the plane of a 2D 
electron gas and the Rashba potential perpendicular to 
this plane. The single particle electron spectrum of thus 
effectively created ballistic quantum (quasi)-lD wire re- 
veals quite a strong deviation from parabolic shape. The 
extent of this deviation depends significantly on the re- 
lation between the strength of confining potential and 
the Rashba SO coupling constant. In order to capture 
the physics of the effective ID wire with nonparabolic, 
spin-splitted spectrum it has been recently suggested to 
consider a modified TL model characterized by two dif- 
ferent Fermi velocities^ vi ^ V2- Without loss of gener- 
ality, this model may be regarded as a multicomponent 
TL liquid introduced by Penc and Solyom^. 

The goal of this paper is to study the persistent cur- 
rents in the class of the multicomponent TL liquids de- 
scribed above, emphasizing the features which can be po- 
tentially observed in experiment. In Sec.^we review the 
derivation of the dispersion relations in the ring geometry 
performed in Refs. .34.43-44 . presenting it in an optimized 
form. In particular, we consider other profiles of the con- 
fining potential which allow to treat the radial part of the 
wavefuction without any further approximation. We also 
estimate the parameter rj — [vi — f2)/(wi -I- W2), which 
measures the asymmetry (nonparabolicity) of the single 
particle spectrum for the recently fabricated semiconduc- 
tor ring with Rashba SO coupling^. In Sec. IIIII we out- 
line the basic steps of bosonization procedure and thus 
introduce our notations. We also quote the main results 
concerning the multicomponent TL liquids. We appre- 
ciate the importance of the zero-mode contributions to 
a persistent current and therefore take a special care of 
the zero-mode part of the bosonized TL Hamiltonian. 
In Sees . Irvl and Ivl we consider persistent currents in the 
cases of canonical and grand canonical ensembles, respec- 
tively. The important issue in these sections is the par- 
ity dependence of persistent currents, which in bosoniza- 
tion picture stems from the topological constraints (selec- 
tion rules) imposed onto the topological excitations (zero 
modes) . We observe the novel features of persistent cur- 
rents caused by the nonzero value of rj. The discussion 
of the results obtained is presented in Sec. IVII 



II. MESOSCOPIC RINGS WITH RASHBA 
INTERACTION: DISPERSION RELATIONS 

The Hamiltonian for a quasi-lD ballistic ring with the 
Rashba SO coupling in polar coordinate sSMM i is a sum 
of the radial, angular and spin-orbit coupling terms 



m 



Tjang 
-"0 



f>^(d^ ld_ 
2m \dr'^ r dr 

2mr^ \ dip 



y(r), (1) 

(2) 



H- ar ( . d 

r \ OLp 



d_ 
'(9r ' 



«aRO-i^^, (3) 



where (7$ — $/$0j 'ti is an effective (band) mass of 
electron, q;r is the Rashba coupling constant, Ur = 
cos ipui + sin Lpa2 , a^p = — sin ipoi + cos ipu2 , and the 
Zeeman interaction is neglected. The confining potential 
V{r) can be modeled by a singular oscillator potential— 
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or by a hard-wall potential 
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In these expressions a can b e asso ciated with a mean 
radius of the ring, and l^ = ^Jh/rmjj or d - with its mean 
width. 

We introduce dimensionless parameters A = a'^/l'^, 
y = r/l^, ao = aB.ma/h^, h = H/huj, and rewrite the in- 
tegrable /i™'' and nonintegrable parts h^^^ = /Jq"^ + ^'^ 
in the new notations: 
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where q = —id^ — (7$. 

We can remove the dependence on (/? in Q (which 
enters through ar and cr<^) by a gauge transformation 

and therefore obtain 
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Note that due to the transformation (jSJ the periodic 
boundary conditions (at zero flux) for an angular wave- 
function have changed to the anti-periodic ones. 



The radial eigenfunctions 
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where L^ is the generahzed Laguerre polynomial, and 



Rn{y) = J-sin\n{n + l){y-VX + l/2) 
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for \/A - 1/2 < y < VA + 1/2, 
elsewhere, (11) 



correspond to the confining potentials Q and (jsj , respec- 
tively. In both cases they are labelled by n = 0, 1, . . . and 
normalized according to 



yRjn{y)Rn{v)dy = (5„i„. 



(12) 



We would like to emphasize that in our consideration we 
do not make use of the approximation A ^ 1, which is 
usually employed when the confining potential is mod- 
eled by a regular harmonic oscillato r'^^i"^^ . This becomes 
an important issue when the mean radius of the ring is 
comparable to its mean width. 
We can find the spectra of JB)) 
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for a singular oscillator and a hard- wall potentials, re- 
spectively, and the matrix elements /i^„ of 10 in the 
basis of Rn{y)' 
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It is remarkable^ that the diagonal elements c„„ identi- 
cally vanish because they correspond to the momentum 
operator in polar coordinates. 

Ideally, we have to diagonalize /i^„ in the infinite basis 
N ^ oo {0 < m,n < N — 1). Practically, it can be only 
done approximately at some large but finite value of N. 
After having diagonalized /ij„n> ^^ can perform another 
- inverse of IJSJ) - gauge transformation 



h: 



diag 



u" — qu' c-1 
"'diag '-^"'diag'-^ i 



5 ^ g-«V'T3/2 ^ ;LAr (18) 



in order to restore the periodic boundary conditions for 
an angular wavefunction. The formal expression (|18|l is 



equivalent to saying that the spin-up e'^^{q) and spin- 
down s'ni^o) dispersion curves in every band n must be 
shifted in 5 by —1/2 and -1-1/2, respectively. 

To describe the main qualitative effects of spin-orbit 
coupling it is sufficient to consider the one-band {N = 1) 
and two-band {N — 2) approximations only. We calcu- 
late the matrix elements (|15|) - (|17ll in the singular oscilla- 
tor basis (|10|l . as well as in the hard- wall potential basis 
(|ll|l . In the former case we have 
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Respectively, in the latter case we have cqi = —cw = 
-8/3 and 
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where fc = 1,2 correspond to a and b, respectively. It is 
easy to compute the integrals in (|24l) and H25|l numeri- 
cally. 

In the one-band approximation we find the dispersion 
relations 

hoo=S±(q) = —{q-q,S,TqRf, (26) 

where f/R = M ^1 -I- 4Q;QaQQ/A6§Q — 1 j is the Rashba an- 
gular momentum. The upper sign corresponds to the 
spin-up branch and the lower sign - to the spin-down 
branch. In this approximation the effect of spin-orbit 
coupling shows up in the relative horizontal shift by 2gR 
of the two parabolic dispersion curves. 

In the two-band approximation the analytical results 
are also available: we can find the solutions of the forth- 
order characteristic equation for 4 x 4-matrix H14|l and 
then shift them according to (|18l) . However, it is easy 
to find the dispersion relations in this case numerically. 
We use the following estimates for the parameters of 
the ballistic ring which has been studied in the recent 
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FIG. 1: The lowest radial band in the two-band approxima- 
tion (qo ~ 3.0). Solid thick line: hard- wall potential. Dashed 
thick line: singular oscillator potential. 



FIG. 3: The lowest radial band in the singular oscillator po- 
tential (qo = 3.0). Dashed thick line: two-band approxima- 
tion. Dashed thin line: one-band approximation. 





FIG. 2: The lowest radial band in the hard-wall potential 
(ao = 3.0). Solid thick line: two-band approximation. Solid 
thin line: one-band approximation. 



FIG. 4: The lowest radial band in the hard-wall potential 
{ao = 6.0). Solid thick line: two-band approximation. Soltd 
thin line: one-band approximation. 



experimenli^: a = 350nni, d = ISOnin, mean free path 
~ 1/xm, m ~ 0.042?7ie. The Rashba couphng con- 
stant OR for the heterostructure used in this experiment 
was « 0.8 • 10~^^eVm. In other heterostructures it can 
reach even larger values (see, e.g. Ref. |47|) in the range 
aR « 1 -^ 6 • lO^^^eVm. With these parameters we can 
find A — 3.78 and ao in the range 2 ~ 12, as well as the 
matrix elements 



aoo = 0.50(0.52), an 
boo = 0.26(0.27), fell = 



= 0.47(0.52), aoi 
0.26(0.28), &01 = 



= 0.11(0.05), 
0.12(0.05), 



and coi = — 0.97(— 2.67) in the singular oscillator (hard- 
wall) potential basis. The results of the band structure 
calculations are presented in Figs. 1 1141 In each figure 
the level e™''(n = 0) is set to zero. We plot the lowest 
radial band SO-splitted into the two subbands (spin-up 
and spin-down) . In Fig. ^ we compare the results of the 



two-band approximation in the hard-wall and the sin- 
gular oscillator potentials for the value of the Rashba 
constant ao = 3.0. In Fig.|21the results of the two-band 
and the one-band approximations in the hard-wall po- 
tential are presented at ao = 3.0. In Fig. O the results 
of the two-band and the one-band approximations in the 
singular oscillator potential are presented at ao — 3.0. 
In Fig. ^ we compare the results of the two-band and 
the one-band approximations in the hard-wall potential 
at ao = 6.0. So, we can see that the dispersion rela- 
tions calculated in the ring geometry appear to be sim- 
ilar to the electron spectra in the wire geometrji22iiMi. 
The common qualitative feature of the wire and the ring 
dispersion curves beyond the one-band approximation is 
their bending, i.e the deviation from the parabolic shape. 
It emerges already in the two-band approximation, and 
including higher bands into consideration would not qual- 



itatively change the picture. The immediate consequence 
of the spectrum bending is the difference in Fermi veloc- 
ities 
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However, due to the time-reversal symmetry expressed in 
the form e^{q) — ei{—q) (at zero flux) we still have 
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Therefore, we have two distinct (in the absolute value) 



Fermi velocities vi 



right 

"ft 



and V2 



VpF . We intro- 
duce the notations vq = (vi -t-f2)/2, Sv = {vi — V2)/2 and 
r] = Sv/vq. For ao — 6.0 we make an estimate rj ^ 0.27 
in the hard-wall confining potential at the Fermi energy 
~ 1.48 measured from e™''(n — 0). 

Recently, there have been studied the effects of non- 
zero 1] on two-point^^ and four-point^ correlation func- 
tions of the TL liquid. In the next sections we would like 
to consider how the spectrum bending discussed above 
modifies the persistent currents. 



III. BOSONIZATION OF THE 

MULTICOMPONENT TOMONAGA-LUTTINGER 

LIQUID 

In this section we present the basic steps of bosoniza- 
tion procedure and the main issues of the TL liquid the- 
ory. In what follows we will use the notations 



27rg Lip 

k = and x = — 

L 2iT 

instead of q and ip in order to employ the standard 
bosonization formulas and to have a transparent relation 
between the description of the ring and the wire geome- 
tries. 

Linearization of the spectrum near the Fermi energy 
yields the second-quantized Hamiltonian 
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where iI)rs{x) and ipLs{x) are the right and the left com- 
ponents of the fermionic field 
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The standard bosonization Ansatz reada^ 
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where 77 = _R(+),L(— ), N^s are the topological excita- 
tions (zero modes) and i^,,s are the Klein factors. One 
can employ the mode expansion of the bosonic fields 
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in terms of bosonic mode excitations b\ {hk^js), a — > 
being a small cut-off parameter. We have also introduced 
the effects of a flixx and a relative shift of the spin-up and 
spin-down dispersion curves by including the boundary 
term 



kj-js — 



fc2 
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into the bosonization formula H31(l . Here ki and ^2 are 
the linearization points for the right spin- up and the right 
spin-down branches, respectively. 

Bosonizing the Hamiltonian (|29() . we obtain 
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The next step is to take into account the standard 
g4c(s) and g2c{s) interactions (forward scattering) which 
specify the Tomonaga-Luttinger modelS4. In our case 
we have the two-component TL liquid, and its continu- 
ous (fc ^ 0) part can be diagonalized by the canonical 
transformation"*^ A = {Aij) 

ihRT,b,^^^,bkRi,bf,^^) = A{dki+,df,^_^_,dki-,di^^_) , 

where A^j = Ay (5'4c(s),52c(s); f^^)- Thus, we obtain 
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where v± is expressed^ through the charge (spin) velo- 
city v^ts): charge (spin) stiffness K^u\ and 6v. 

We note that the zero-mode part in (|35|l (the second 
and the third lines) cannot be fully diagonalized by the 
continuous canonical transformation A, since the charge 
and spin current excitations 

Jc = NR^-~N^-+NR^~NLi, (36) 

Js = Nr^-^Nl^^-Nr^+Nli, (37) 

as well as the charge and spin number excitations 

Nc - Nr^ + NL^ + Nri + Nli , (38) 

Ns = Nr^ + Nli; - Nr^ - Nli (39) 

take integer values, and therefore they can be only trans- 
formed by a modular transformation which maps a grid 
onto itself. For this reason the grand partition function 
So corresponding to zero modes cannot be expressed in 



terms of the standard Jacobian theta functions^^ as it 
occurs in the spinless caseSZi and even in the spinful case 
with spin-orbit coupHngS. It is exactly the spectrum 
asymmetry parameter Sv that violates the factorization 
of So into different Jc, Js,Nc, Ng sectors. In the canonical 
ensemble the factorization of the partition function Zo is 
still possible, but the terms proportional to 5v would also 
cause some corrections. 

In 135|l we have also introduced 



Jc = Jc- 4:q<s>, Js ^ Js - 4(5q, 
Nc = Nc- 4go, N, = iV„ 
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The similarity sign turns into the equality sign for the 
one-band approximation, i.e. when we deal with the 
parabolic dispersion relations. The first expression in 
(|40|l also shows that the entire dependence on the mag- 
netic flux is contained in the zero- mode part of the Hamil- 
tonian H35|l . and therefore it is sufficient to calculate 
ZqC^q) in order to find the persistent currents in the 
(grand) canonical ensemble. 



IV. PERSISTENT CURRENTS: FIXED 
NUMBER OF ELECTRONS 



The curly brackets in {Jc,Js} mean that the certain 
topological constraints are imposed onto the pair of val- 
ues Jc, Js- These constraints, or selection rules, are 
different for different parities of the electron number. 
This gives rise to the parity dependence of persistent 
currents^-"'. Before discussing it in further details, we are 
already able to state that the whole effect of the spec- 
trum asymmetry in the case of fixed number of electrons 
is the renormalization of the Rashba angular momentum 

The topological constraints in question have been de- 
rived in Ref. |23 for spinless TL liquid, and the case of 
the spin-i TL liquid has been considered later on by dif- 
ferent authors^LSiiS. In what follows we will use the 
formulation of Ref. 152. 



A. iVe = 47Vo + 2 and iVe = 47Vo 

When the number of electrons in the ring equals to 
Nc — 4A'o -I- 2, we have the following possible combina- 
tions of {Jc, Js}'- 



{AucAns}, {4nc + 2,4n, + 2}. 



(47) 



The summations over n^ and Ug are unconstrained and 
run from — cx3 to cxd. Using the properties of the Jaco- 
bian theta functions (Appendix^, we can find that the 
persistent current /'^'(g^) equals to 



In case of fixed number of electrons persistent current 
equals 



/($) = -^A^(*), 



where AF(<I>) is an oscillating part of the free energy. 

In the following we adopt the units in which $0 = 1 
and introduce the characteristic scale for a current Jo = 
47rvo/i. 

In the canonical ensemble the charge and spin number 
excitations are forbidden: Nc = Ng = 0. Therefore, 
different topological sectors disentangle, and we find 
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summing over {J^ Jg} sectors. We have introduced the 
Hamiltonian 
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(45) 
along with the effective Rashba angular momentum [cf. 

m] 






Sq 



2go-l 
2vgKg 



Sv. 



(46) 



ttT' 



y3{T^q<i>,ic)d3{'^qK'^ ,1s) + d'4{Trq<s>,ic)04{T^qn'^ ,1s) 



h{T^q'i',lc)0?,{T^q^/ ,ls) + 0i{nq^,-ic)9i{'Kq^^ ,-lgy 



(43) where 



lc(s) = e 



-ttLT lv^(^)K^ 



(48) 



(49) 



In the case 5u = we find the agreement between (|48|) 
and the respective expression of Ref. |33. Considering 
further the noninteracting limit 



7c = 7, =7 = e-^^^/''«. 



(50) 



we can simplify (|48|) down to 

nw N v^ 47rT(-I)" , ^ , 

I^ H9* = 2^ ■ u/o 2 rr,r ~, cos{2TmqYi)sm{2Trnq^), 
^-^ smh ZTT^nJ /In) 
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_ (51) 

and recover the result of Ref. |33. 

For Ne = 4:Nq we can find the persistent current by 
the mere shift of the argument 



/(2)(q$)=/W(q$ + f/2), 



(52) 



which is valid in both noninteracting^^ and interacting^ 
considerations. 



B. Ne = 4Afo + 1 and Ne = 4iVo + 3 

When Ne = 47Vo + 1 or iVg = 4iVo + 3, i.e. for the 
odd electron number, we have to sum in (|44|l over the 
foUowing combinations of {Jc, J,}: 

{An^ + 1, An, + 1}, {4nc + 1, An, + 3}, 

{Anc + 3, An, + 1}, {An^ + 3, An, + 3}. (53) 



Due to 5v ^ Q different topological J and N sectors re- 
main entangled. We can rewrite \&Hq{<^ , fi) / Iq (up to 
some irrelevant additive constant term) in the form 

Ac(Jc - z,^f + 2ri{Jc - z^){N, - zb) + ^siN, - z^f 
+ K{Js - ZR,)2 + 2ri{J, - ZR)(iVc - z^) + Vc{N^ - z^,f 

with 77 and /q introduced in Sees. HTl and II VI respectively; 



We can find then 






47rr 



■ sin(47rng$). (54) 



■^ sinh(7rnLT/wc-?^c 

It is remarkable that the result does not depend on q^ , 
and therefore the spin-orbit coupling does not show up 
in this case at all. In the noninteracting limit v^ = vn 
and Kc = 1, and we are again in agreement with Ref. l33l 
From this perspective the formula (13) of Ref. 1371 seems 
to be incorrect. The probable source of error might be in 
the following: even though all possible topological con- 
figurations of the ground state were correctly established 
(equivalently to H53|) '). not all of them were included in the 
computation of the free energy. Therefore, the statement 
of Ref. I33 that for finite temperatures there might occur 
the spontaneous currents (at zero flux) is misleading. In 
contrast, according to (|54|l we always have /'■^■'(O) = 0. 



V. PERSISTENT CURRENTS: FIXED 
CHEMICAL POTENTIAL 

Let us now consider the case when the mesoscopic ring 
is weakly coupled to the reservoir with chemical potential 
11. In this situation the charge and spin number excita- 
tions are allowed. For this reason the terms mixing J 
and N sectors in (|35|l do not vanish, and therefore the 
consequences of the spectrum asymmetry {5v 7^ 0) may 
become more diverse. 

In order to find the persistent current 



^('i>,Ai) = -^Afi($,/x), 



(55) 



we need to calculate the oscillating part of the thermo- 
dynamic potential 



Af]($,/i) = -rin 



E 



^ffo(*,M)/T 



(56) 



{.Jo,Nc,,Js,N,} 



where the curly brackets in {Jc^Nc^Js^Ns} again mean 
the topological constraints, and 



Ac(s) = 



Vc(s)Kc(s) 



'Ms) 



Vc{s) 



K 



{s)VO 



and 



z$ 



Aqi 



ZR^AiSq- 



ZB = 0, 



= 4Uo 



??(./m + ^c) 
2(^cA, - r;2) 

(i^cAs - ry2) 



(58) 

(59) 
(60) 



1 measures the 



In the above expressions /^ — ^^ 
difference between the chemical potential and the Fermi 
energy, and zb is some argument, which would have been 
non-zero, had we included the Zeeman interaction. 

We remark that the charging energy Ec and the gate 
voltage Vg are not explicitly included in ifSTI) . We infer 
that they can be effectively taken into account by redefin- 
ing Ac ^ Ac -I- 16Ec/Iq and fi ^ ^ + eVg, respectively. 

In Appendix ^ we introduce auxiliary functions 
G'(2$, Zb, zr, z^) and G"(z$, zb, zr, z^^), which allow to 
conveniently express S(<I>, /i) = g^-^^{'^^t^)/T and /($,/j,) 
for different electron numbers N^ in the ground state (at 
T = and <I> = 0). In the noninteracting limit and/or 
in the limit 77 = the functions G and G" can be sim- 
plified and rewritten in terms of the standard Jacobian 
theta functions (due to the relations quoted in Appendix 
IX|) . In general, they are expressed in terms of the Siegel 
theta functions^^ (i.e. theta functions of higher degree), 
and therefore it is not straightforward to establish their 
asymptotics. However, for small enough temperatures 
the functions G and G" can be easily computed numeri- 
cally. 



A. Ne= ANfj + 2 and N^ = ANq in the ground state 

The selection rules formulated in Ref. |52 are also ap- 
plicable in our consideration. Thus, for ATe = ANq + 2 
electrons in the ground state we have to choose those J's 
and iV's which satisfy the following requirements: 

1) Jc and Nc are either simultaneously even or simul- 
taneously odd. The same holds for J, and Ns\ 



Hoi^.lA 



AL 



V^KcJr + V.KsJ^ 






+ 2Sv ( JcNs + JsNc + AvqNc + ASvJ, 



M^c (57) 



2) Jc±Js + Nc± Ns takes values . . . , -4, 0, 4, . . . 

In fact, this concisely formulated prescription as- 
sumes the summation over 16 different combinations of 



{Jc,Nc,Js,Ns} (see Appendix IB|) . Having properly per- 
formed it, we obtain 



and 






/W($,M) = 27rT=i('^'^^ 



Si($,m)- 



(61) 
(62) 



(63) 



In the noninteracting limit Xc(s) — Vc{s) = 1 the per- 
sistent current (|63(l equals to 



47rT^(-l)^ 



sin ^ z^ cos ^ z'^ sin ^z'^cos^z'^ 



where 



sinh 



2$ 
Z<[) 



Io{l+v) 



ZB, 
ZB, 



sinh 



= ZR + Zf, 
= ZK - Za 



Io{l-v) 



(64) 

(65) 
(66) 



Recalling that zb = 0, we can easily establish the Fourier 
coefficients /„ of the series 



(|M|l = 2_, ^n sin 27rng$ . 



(67) 



In the limit T — > we find that /„ equals to 



2/r 



(-1)" 



nn 



nn 



nn 



cos -— ZR cos -—Zf^ + r] sm -— zr sm -— z^ 



So, we obtain additional modulation due to 77 7^ 0. 

In the limit rj = the grand partition function l|61|) is 
proportional to 



SMf 



z<s,,aA 9, 



i=l 



with 



9,(|zR,«^)e,(|z„/3,4)0,(|zB,/3,4) 

(68) 



a, 



c(s) 



JO^c(s) 



Pc{s) 



JO-'c(s) 



(69) 



The coefficient of proportionality in 1)68(1 is omitted, since 
it does not depend on z$ and therefore becomes unim- 
portant for the calculation of the persistent current H63|) 
in this limit. We note that the expression H68|l coincides 
with that derived in Ref. |33 

In Figs. 03 the persistent currents for the ring with 
Ne = 4A'o-f 2 electrons in the ground state (i.e. at T = 0) 
are presented at T = O.OOS/q. To parameterize the e-e in- 
teractions we use the relation between the parameters of 
the TL model and the Hubbard model at low densitie^: 
Ac = As = 1, I'c = 1 + u, j^s = 1 — u, where u — U/ttvq 
and U is an interaction parameter in the Hubbard model. 
Other parameters are 6q = 0.15, go = 0.04, /^ = 0.11. 

In Figs. Inland El we show how the persistent current in 
rings with interacting (u = 0.4) and noninteracting (u = 
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FIG. 5: Persistent currents I/Io vs. flux $/$o at u — 0.4 and 
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FIG. 6: Persistent currents I/Io vs. flux ^/^o at m = 0.0 and 
different values of rj (fixed p, Ne = 4:No + 2). 
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FIG. 7: Persistent currents I/Io vs. flux $/<l?o at r; = 0.3 and 
different values of u (fixed fi, N^ = 4A'o -I- 2) . 



0.0) electrons is modified when we vary the spectrum 
asymmetry parameter 77. One can see that in the first 
case increasing 77 develops a paramagnetic "tooth" near 
g$ = 0, while in the second case altering rj changes the 
size of the diamagnetic "tooth" near q$ = 1/2. 

The dependence of the persistent current at 77 = 0.3 
on the interaction parameter u is depicted in Fig. \7\ 

Of course, an alteration of other parameters Sq, qq, 
fn may also change the picture of the current, but it is 
obvious that the modification due to 77 is remarkable. 

The results for Nf. = 47Vo are obtained by the shift in 
either z$ or z^ (see, e.g. Ref. Isij) : 



^2(*, m) = G{z^ + 2, ZB, z-R, Zi_,) 
= G{z^,ZB,z-R,Zf, + 2), 



(70) 



and similarly for /'^^^(<f>, /i). 




FIG. 8: Persistent currents I/Io vs. flux $/<l?o &t u — 0.4 and 
different values of rj (fixed /i, odd N^). 



B. Ne ^ 4No + 1 and Ne = 4Afo + 3 in the ground state 

First of all, wc would like to remark that the ground 
state with the odd electron number {Ne = 4:No + 1 or 
A^e = 4A^o + 3) can be justified only for repulsive in- 
teractions at very low temperatures: T <C U/L ~ Iqu 
(see, e.g., Ref. |52for the discussion). In this case we can 
formulate the selection rules in terms of the function G 
and thus obtain the following expression for the grand 
partition function: 

S3($,/i) = G{z^ + l,ZB,z^, + l,Zi_,) 

+ G(Z<I. + 1,ZB,ZR - 1,Z,,): (71) 

%{^,^) = G"(z<i. + 1,ZB,ZR + l,Zp) 

+ G'{z.s, + l,ZB,ZK-l,Zf,). (72) 



Respectively, 



/(^)($,^) = 27rT 



^3(*,/")' 



(73) 



In the noninteracting limit the grand partition function 
(|7T|l is proportional to 



where 7± = e wJ^vJ and 



„ _ Zg, + ZR ZB + Zf, „ 
Zfh R — ::: i :: , z-c 



^$,R 



"' B,/^ 



Z<S> - Zr Zb - Z^ 



In the limit 7; = the grand partition function H71f) is 
proportional to 

Oi (2 ^*'"c) 6'3 \^ZB,Pfj 6*4 (2 ^R '"^) ^3 (2 ^A''/^; 
+ 6I3 (^-z$,a^j 6*4 (2 ^B,/3|) ^3 (2 ^R'"^) ^4 (2^'^''^< 



If we put in this expression zr = zb = (no spin-orbit 
coupling and no magnetic field) and as = (3s, we obtain 
the formula equivalent to that derived in Ref. ISlt 



S3($,/i) 



73 [ -;zz<s>,ac } "74 



(f 



4 
z$, a^ 






(74) 



In Fig. IHlwe present the persistent currents at ti = 0.4 
and different values oirj, other parameters being the same 
as in the even-iVg case. One can observe that in the pres- 
ence of r; ^ the persistent current is modified consid- 
erably. In particular, the effects of SO coupling show 
up in the grand canonical ensemble in the ring with odd 
number of electrons in the ground state, in contrast to 
the ring with fixed odd number of electrons where such 
effects are absent [cf. (|54|) ]. Due to 77 ^ there also 
happen transitions between the diamagnetic and param- 
agnetic behavior at g$ — 0, ±1/4, ±1/2, . . ., and for some 
particular value of r] (which is close to 0.45 for u — 0.4) 
a period-halving occcurs (i.e. the period becomes equal 
to q^/A). 



VI. CONCLUSIONS 

We have demonstrated that the account of the higher 
radial subbands leads to the modification of the spectrum 
of electrons in the ring geometry. In particular, the dis- 
persion relations obtained are characterized by two dif- 
ferent Fermi velocities. Therefore, in order to describe 
the joint effect of the e-e interactions and spin-orbit cou- 
pling in ballistic rings it is reasonable to consider the 
multicomponent Tomonaga-Luttinger model. The inter- 
play between the spin-orbit splitting, nonparabolicity of 
the spectrum (rj y^ 0), e-e interactions and parity effects 
is reflected in the spectral properties of persistent cur- 
rents. We have studied this within both canonical and 
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grand canonical ensembles. Below we summarize the ba- 
sic features of persistent currents caused by the subband 
nonparabolicity. 

In Sec. IIVI we have shown that in the rings with the 
fixed electron number persistent currents for 77 ^ do 
not change a lot compared to the case rj — 0. The reason 
for that is the specific selection rule dictated by the par- 
ticle conservation in the system. On the other hand, in 
the rings with fixed chemical potential we have observed 
considerable modifications. Thus, for the even number of 
electrons in the ground state we have obtained the gen- 
eration of new harmonics and complication of the cur- 
rent shape. We have observed how the picture changes 
if we vary the effective Coulomb interaction parameter 
u, while fixing the effective nonparabolicity parameter 
r/, and vice versa, and found these changes remarkable. 
The modification of the current in the case of odd num- 
ber of particles in the ground state appears to be even 
more drastic: varying rj we can perform the transition 
from diamagnetic to paramagnetic behavior at the flux 
values $ = 0, ±$o/4, ±$o/2, . . . This makes the current 
essentially different from that derived from the parabolic 
dispersion relations for the same parity. 

We therefore conclude that the nonparabolicity of the 
single particle electron spectrum produces a deep entan- 
glement of the e-e repulsion and spin-orbit coupling pa- 
rameters, and this leads to the new features of the cur- 
rent, which can be detected experimentally. 

Recently, there has been the certain development in 
fabrication of toroidal carbon nanotubes^, and in the 
study of persistent currents in such systems^^. Since the 
electronic spectrum in carbon nanotubes deviates consid- 
erably from the parabolic shape, we expect our results to 
be applicable in such realization as well, provided the 
selection rules are properly modified. 
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APPENDIX A: PROPERTIES OF THE 
JACOBIAN THETA FUNCTIONS 



We use the following definition of the Jacobian theta 
functionsSii 



6)2(2,7) = 27i/4^7"("+i)sin(2n+l)z, (A2) 



n=0 



93(2,7) = l + 2^7"'cos2nz. 



(A3) 



ra=l 



614(2,7) = l + 2^(-l)"7"'cos27iz. (A4) 

n=l 

Note that there exists the alternative definition which 
uses the different argument: 2: — > ttz. 

Functions ^3 and 64 are periodic under the shift z -^ 
z + TT, while 01 and 02 change their signs. Note also that 

^3(2+1,7) =Oi{z,-f), ^2(2+1,7) =-9i{z,-f). 

Making Poisson resummation, it is easy to prove the 
useful formulas 



E^-' 



{k+zY _ 



k— — oo 



J2 {-!)'' e-^"'-'' 
fc— — 00 

From the relations 



+^r - 



^9s{7Tz,e-^'^^),{^5) 



/!^^(^ 



-TT^'/a 



. (A6) 



63.2 (2z, 7^) = - (03(^, 7) ± Oi{z, 7)) , (A7) 
it is easy to deduce that 

63 (221,7^) 03 (222,72) + O2 (221,7^) 02 (222,72) 
== 2^3(^i,7i)6'3(z2,72)+ 2^4(21,71)^4(^2,72); (A8) 

03 (221,74) 03 (222, 7I) - O2 (221, 7*) 02 (222,72') 

= 2^3(^1' ■^1)^4(2:2, 72) + 2^4(21,71)^3(22,72)- (A9) 
Moreover, for 71 = 72 = 7 

[1 = 03 (21 + 22, 7') ^3 (21 -22, 7'), (AlO) 
[3 = 04 (21 -h 22, 7') 04 (21 -22, 7')- (All) 

There also holds the following identity 

4 
^0i(zi,'y)0i(z2,'y)0i{z3,j)0,{z4,j) 

i=l 

= 203(2^,7)^3(2^', 7)^3(2^,7)03(2^', 7), (A12) 
where 



„ _ 21 +22 , 23 + 24 „ 

2i 9 — :: ^ ::: , 2q 



'1,2 



"' ^3,4 



21 - 22 23 - 24 



0i(z,7) = 27i/4^(-l)"7"("+i)sin(2n+l)2,(Al) 



n=0 



In Ref. I53 one can also find useful expressions for the 
logarithmic derivatives 0^(2,7)70^(2,7). 
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APPENDIX B: TOPOLOGICAL CONSTRAINTS 
FOR FIXED CHEMICAL POTENTIAL 

The topological constraints formulated in Ref. ^2 for 
Ne = 4iVo + 2 electrons in the ground state lead to the 
following possible combinations oi {Jc,Nc, Js,Ns}: 

{4nc,4mc,4ns,4ms}, 

{4nc + 2, Anic + 2, An^ + 2, 4ms + 2}, 

{4nc + 2, 4mc + 2, 4n„ 4m J, 

{4nc, 4mc, 4^^ + 2, 4m, + 2}, 

{4nc, 4toc + 2, 471,, 4to, + 2}, 

{4nc, 4toc + 2, Aus + 2, 4to4, 

{4nc + 2, 4toc, 4ns, 4™, + 2}, 

{4nc + 2, 4toc, 4ns + 2, 4™^}, 

{4nc + 1, 4toc + 1, 4ns + 1, 4ms + 1}, 

{4nc + 3, 4toc + 3, 4ns + 3, 4ms + 3}, 

{4nc + 1, 4toc + 1, 4ns + 3, 4ms + 3}, 

{4nc + 3, 4toc + 3, 4ns + 1, 4ms + 1}, 

{4nc + 1, 4mc + 3, 4ns + 1, 4ms + 3}, 

{4nc + 1, 4mc + 3, 4ns + 3, 4ms + 1}, 

{4nc + 3, 4toc + 1, 4ns + 3, 4ms + 1}, 

{4nc + 3,4mc + l,4ns + l,4ms + 3}. (Bl) 

Therefore, we have to perforin 16 different unconstrained 
summations over Uc^ mc, Us, mg from — oo to cx). 

Let us define for convenience the functions [cf. Ref. \5^ 



/i(z$,2;b) 



f2izR,Zfj,) 



E --'" 



-Iohi{z^,ZB)/T 



TlcTTlg — — oo 



E -^'"' 



2(zn,z^)/T 



(B2) 



(B3) 



Us ,mc — — oc 



where 



z$ 



ZB 



/ii(z$,zb) = Xcinc-—] +iys[rns~ — 
+ 2rj (n, - ^) (ms - ^) , 

h2{zR, z^) = As (us - -^j + Vc [nic - -^j 
+ 2ry {ns - ^) (m, - ^) , 

and we must demand Vs^c — ff' > and Vc\s — 77^ > to 
ensure the convergence of the series in (jB2ll and (|B3ll . 
Note that the functions /i and /2 have period 4 in each 
argument. 

If we introduce 

g(z$,ZB,ZR,Zp) = [/i(z$,zb) + /i(2$ + 2,zb + 2)] 
X [h{zK,z^) + f2{zR + 2,z^ + 2)] 

+ [/i(2:<i>,ZB + 2) + /i((z<i, + 2,ZB)] 
X [/2(2:r,Zp + 2) + /2(zr + 2,Z^,)], 



then it becomes easy to see that the grand partition func- 
tion equals to 

G{z^,zb,zr, Zfi) = g{z<^, ZB, ZR, Zf^) 

+ g{zg. + l,ZB + l,zn + l,z^, + l). (B4) 

We can also define 
f[{z>^,ZB) = ^ E [^c (4nc - z$) + ?7 (4ms - zb)] 



xe 



-Iohi(zis,.z-B)/T 



(B5) 



as well as g' and G' replacing /i by /{ in the above defi- 
nitions. 

In the nonintcracting limit 



, , 1 + ?? / , Z^+ZB 

hi{z^,ZB) = —^ — (nc-t-m. 



2 

1-7? 



Z<S> - ZB 



So, we can transform Uc, nis by a modular transforma- 
tion. Similarly, we can proceed with h2 and Us, mc- 

In the hmit r] ~ the sum /i(z$, ZB)-f/i(z$ + 2, zb + 2) 
is proportional to 

6*3 ijZ^, ac) 03 f -ZB, /3s j + 04 ( T^*> "cj 04 f -2b, ^sj , 

while /i(z$ + 2, zb) + /i(z$ + 2, zb) is proportional to 
6*3 (^-z^,acj 04 (^ZB,/3sj +03 (^z*,acj 04 (^zb,^sJ , 

where ac and /3s are given by l|()9|) . Similar relations 
hold for /2. Using ljA8|l and (|A9|) . one can establish the 
expression (|68|l . 

For the odd number of electrons in the ground state 
one can also establish the selection rules and calculate 
the grand partition function. In terms of the function 
G it is presented in (|71|l . An alternative way to express 
S3(4>, ^) is to introduce 

G{z^, ZB, zr, Zp) = g{zq,, ZB, zr, z^,) 

+ g(z$ + 1,zb + 1,zr + l,Zp + 1), (B6) 



where 



g(z$,ZB,ZR,2;^,) = .g(z$ -I- 1,zb,zr-|- l,z^) 
+ g{zq, + 1,zb,zr- 1,z^) 
= /i(z$,zb)/2(zr,z,,) (B7) 



and 



/i(z$,zb) == /i(z$ + 1,zb) + /i(z$ + 1,ZB + 2) 

+ /i(z$ + 3, Zb) + /i(z<i, +3,ZB + 2), 

/2(ZR,Z,,) = /2(zr + 1,Z^) +/2(zR + l,Z^ + 2) 

+ /2(zR + 3,Z^)+/2(zR + 3,Z^ + 2). 
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One can show that the functions /i and /2 have period 
2 in each argument and equal to 

oo 

fi{z^,ZB) = J2 e-^«''^(^-^-)/^, (B8) 

nc^'^Tls — — oo 
oo 
MZK,Z^) = Y. e-^oS.(.K,..)/T^ (gg) 



Us ,mc — ~OQ 



where 



/ii(z$,zb) 



z^ + I 



j + i^s (r 



2^7 ("c- ^ — j (ms-Y 



ZB 



(BIO) 



h2izB.,z^) 



+ 2rj n. 






In the noninteracting hmit 

I , . l + r] f~ , ~ z^ + ZB + l 

/ii(z$,zbJ = — - — 1 ric + m. 
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rir — rrio — 



z$ - 2;b + 1 



and we transform fic, rhg by a modular transformation. 
Similarly, we proceed with /12 and fig, fa^. 

The limit 7] = for /i and /2 is straightforward. 
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